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Viral infection depends on a complex interplay between host and viral factors. Here, we link host 
susceptibility to viral infection to a network encompassing sulfur metabolism, tRNA modification, 
competitive binding, and programmed ribosomal frameshifting (PRF) . We first demonstrate that the 
iron-sulfur cluster biosynthesis pathway in Escherichia coli exerts a protective effect during lambda 
phage infection, while a tRNA thiolation pathway enhances viral infection. We show that tRNA Lys 
uridine 34 modification inhibits PRF to influence the ratio of lambda phage proteins gpG and gpGT. 
Computational modeling and experiments suggest that the role of the iron-sulfur cluster 
biosynthesis pathway in infection is indirect, via competitive binding of the shared sulfur donor 
IscS. Based on the universality of many key components of this network, in both the host and the 
virus, we anticipate that these findings may have broad relevance to understanding other infections, 
including viral infection of humans. 
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Introduction 

All viruses require host translational machinery in order to 
replicate. To make use of this machinery while avoiding 
detection, viruses have evolved a number of unconventional 
translational strategies (Gale et al, 2000), including internal 
ribosome entry sites (Jang et al, 1988), leaky scanning 
(Schwartz et al, 1990), ribosome shunting (Schmidt-Puchta 
et al, 1997), reinitiation (Hemmings-Mieszczak et al, 1997), 
and programmed ribosomal frameshifting (PRF) (Gesteland 
and Atkins, 1996). Many of these strategies are unique to 
viruses, and thus are potential antiviral targets (Harford, 
1995). 

During PRF, the translational machinery slips backward or 
forward by one nucleotide while decoding a so-called 'slippery 
sequence,' a heptanucleotide sequence of the form XXXYYYZ 
in the mRNA transcript (Figure 1). Outside of PRF, frameshift- 
ing is a rare event, occurring fewer than 5 x 10~ s times per 
codon (Kurland, 1992). PRF can increase the occurrence of 
frameshifting by several orders of magnitude, depending on 
the slippery sequence itself as well as on the presence or 
absence of adjacent secondary structures that favor transla- 
tional pausing and increase frameshift frequency (Gesteland 
and Atkins, 1996). In many instances of PRF, slippage causes 
the translational machinery to bypass a stop codon just 
downstream of the slippery sequence, resulting in the 
expression of an elongated protein (Gesteland and Atkins, 
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1996). Therefore, normal expression from these transcripts 
results in two products, a shorter product based on the initial 
translation frame plus a longer product, whose ratio depends 
on the frequency of slippage and is thought to have evolved to 
a value optimized for virus production (Shehu-Xhilaga et al, 
2001;Dulude et al, 2006). 

PRF is an evolutionarily conserved phenomenon observed 
in a range of viruses, from bacteriophage to notable 
human pathogens (Table I; Baranov et al, 2002) . For example, 
in HIV, the Gag to Gag-Pol expression ratio is regulated by 
PRF and is necessary for the proper loading of reverse 
transcriptase into the new viral capsids; disturbing this 
ratio can inhibit viral replication (Shehu-Xhilaga et al, 2001; 
Dulude et al, 2006) . The role of PRF in the SARS life cycle is less 
clear, but is necessary for expression of its RNA-dependent 
RNA polymerase (Plant et al, 2010). In Escherichia coli 
bacteriophage lambda, PRF regulates the expression of 
proteins gpG and gpGT (Levin et al, 1993). As lambda phage 
genes G and The between its major tail gene (V) and its tape 
measure gene (H), they are thought to have a role in tail 
formation (Xu, 2001) . However, proteins gpG and gpGT do not 
make up any of the final tail structure (Levin et al, 1993), and 
thus may act as chaperones or scaffolds for gpH, assisting in 
the assembly of the tail shaft around the tape measure protein 
(Xu, 2001). 

Reading-frame maintenance is a common function of tRNA 
modifications (Urbonavicius et al, 2001). tRNAs contain 
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upwards of 80 known modified nucleosides (Limbach et al, 
1994) , and many of these are strongly conserved across species 
(Sprinzl and Vassilenko, 2005) . However, their loss is generally 
not lethal (Alexandrov et al, 2006; Phizicky and Hopper, 2010) . 
In E. coli, four tRNA modifications require the addition of a 
thiol moiety and are made via two functionally distinct 
mechanisms (Nilsson et al, 2002). The first mechanism 
requires proteins that contain iron-sulfur (Fe-S) clusters and 
includes 2-thiocytidine formation at position 32 and 
methylthio formation of ms2i6A at position 37 (Jager et al, 
2004; Pierrel et al, 2004). The second mechanism is 
independent of Fe-S clusters and includes 4-thiouridine at 
position 8 and 2-thiouridine at position 34 (Ikeuchi et al, 
2006) . For each of these modifications, sulfur is sequestered by 
E. coli's primary cysteine desulfurase, IscS. Sulfur is then used 
for one of these four tRNA modifications or for Fe-S cluster 
biosynthesis. In addition to their roles in tRNA modification, 
Fe-S clusters are critical cofactors involved in many cellular 
processes including respiration, central metabolism, environ- 
mental sensing, RNA modification, DNA repair, and DNA 
replication (Py and Barras, 2010). 

We previously identified two pathways downstream of 
IscS in a screen to identify E. coli genes with a significant effect 
on lambda phage replication (Maynard et al, 2010). Specifi- 
cally, lambda phage replication was inhibited following 
deletion of several members of the 2-thiouridine synthesis 
(TUS) pathway leading to 2-thiouridine modification of 
tRNA Lys/Glu/Gln , a pathway independent of Fe-S cluster 
biosynthesis. Conversely, we found that knocking out mem- 
bers of the Fe-S cluster biosynthesis (ISC) pathway enhanced 
lambda replication. 

The current investigation was therefore motivated by this 
overarching question: how do the host's TUS and ISC 
pathways control viral replication? We were particularly 
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Figure 1 Programmed ribosomal frameshifting. Schematic of a -1 
programmed ribosomal frameshift. P-site tRNA slips in the -1 direction at the 
'slippery sequence.' 
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Table I Viruses dependent on programmed ribosomal frameshifting 



Virus 


Slippery sequence 


Region 


HIV-1 


UUUUUUA 


gag/gag-pol 


SARS-CoV 


UUUAAAC 


ORFla/lb 


Lambda phage 


GGGAAAG 


G/GT 



intrigued by the possibility that these pathways interact. 
Furthermore, the substantial evolutionary conservation of 
both host genes and viral PRF mechanism strongly suggested 
that our observations would have relevance beyond E. coli and 
lambda phage. Our current observations shed light on a 
complex network that extends from host metabolism and 
tRNA modification to viral translational regulation and finally 
to virion production. We find that 2-thiouridine hypomodifica- 
tion of t RNA Lys/Glu/Gln causes increased translational frame- 
shifting, changing the ratio of the critical lambda phage 
proteins gpG and gpGT We also show that IscU is linked to gpG 
and gpGT expression by competitive inhibition of the TUS 
pathway. Knocking out core members of the ISC pathway 
increases lambda phage replication by relieving competition 
for sulfur, allowing the TUS pathway to increase the rates of 2- 
thiouridine formation, thus reducing frameshifting. Targeting 
tRNA modifications to alter frameshifting rates, to which viral 
structural proteins are likely to be uniquely sensitive, presents 
a novel antiviral strategy. 



Results 

Viral replication can be slowed by deletion 

of TUS and accelerated by deletion of ISC genes 

In a previous genome-wide screen (Maynard et al, 2010), we 
found that plates growing TUS pathway deletion strains 
[AtusA, AtusE, and AmnmA) infected with lambda phage 
produced plaques of an unusually small diameter compared 
with wild-type (WT) E. coli. In contrast, several ISC pathway 
deletion strains (AiscU, AhscA, and AhscB) produced 
abnormally large diameter plaques compared with WT. 

We investigated viral replication in these strains more 
thoroughly by culturing each strain in liquid media, in the 
presence and absence of virus. Without phage, WT E. coli grew 
in exponential phase for several hours, slowing and reaching 
stationary phase as available nutrients decreased (Figure 2A, 
blue line). In the presence of lambda phage, the culture 
exhibited three phases (gray line) . First, the culture underwent 
exponential growth similar to the uninfected culture. After 
~ 6 h in the WT strain, viral lysis overtook the culture and the 
culture began to clear (lytic phase). Lambda phage is a 
temperate phage with both lytic and lysogenic life cycles 
(Landy and Ross, 1977); as a result, after a roughly 3-h lytic 
phase, lysogenized bacterial population growth took over the 
culture. 

The time courses of lambda phage infection of the tasA and 
iscU knockout strains (the first members of the TUS and ISC 
pathways, respectively) differed significantly from that of WT. 
The AiscU culture entered the lytic phase very quickly (red 
curve, Figure 2B), while the infected AtusA strain grew in 
exponential phase substantially longer than the WT strain, 
leading to a higher turbidity before the lytic phase occurred 
(green curve). We therefore decided to classify strains that 
entered the lytic phase earlier as more susceptible to viral 
infection (red) and those that entered the lytic phase later as 
less susceptible (green). This susceptibility classification 
was consistent with our earlier plaque assay experiments 
(Maynard et al, 2010). 
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Figure 2 Dynamics of lambda phage infection. (A) A typical trace of a WT £ coli culture infected with lambda phage. Several hours after infection, lysis begins to 
outpace £ coli growth and absorbance begins to decrease. Regrowth is due to the lambda lysogen population. (B) Infection dynamics of infected AiscU and AtusA 
cultures were compared. (C) The displacement from the WT growth curve. In (C), the bars indicate the 95% confidence interval (CI). For (A) and (B), absorbance was 
recorded over the course of 1 6 h for three biological replicates with four technical replicates for each biological replicate (**P< 0.01 , P-values were calculated using an 
independent two-sample f-test). 



We then defined a metric to compare the resistance of 
different strains with infection with phage lambda. One useful 
way to summarize the information contained in a time course 
is to consider each experiment as a vector in n-dimensional 
space, where n is the number of time points taken. In earlier 
work, we found that it was useful to normalize these infection 
time course experiment vectors by the growth rate (Maynard 
et al, 2010). However, we have since found that bacterial 
growth and viral growth are tightly linked, and our attempts to 
numerically correct for or eliminate the effect of bacterial 
growth rate often distorted the raw data. Bacterial growth is 
therefore an important (but not sufficient) consideration in 
classifying our strains. 

As a result, we determined our new metric simply by 
calculating the Euclidean distance between the knockout and 
WT strain time course vectors under infection conditions. 
Combining this distance with our susceptibility classification, 
we obtain a displacement vector that points in a negative 
direction for strains that are less resistant to virus and a 
positive direction for more resistant strains. The displacement 
vector, which we call the 'host resistance', enables us to 
compare the magnitude and direction of each strain relative 
with the WT strain and with each other. For example, the AiscU 
strain is significantly less resistant to viral infection than the 
WT strain, while the AtusA strain is more resistant (Figure 2C) . 



Deletion mutants in Fe-S biosynthesis and tRNA 
thiolation exhibit altered viral infection dynamics 

As both the TUS and ISC pathways are linked to sulfur 
metabolism through IscS, we wondered about the interplay 
between these two pathways and what might be the 
mechanistic basis for their effect on lambda replication. We 
began our investigation by examining the impact of other 
known ISC and TUS pathway members on lambda phage 
replication (Figure 3A; see Supplementary Figure 1 for line 
plots of individual replicates) . First, we investigated the effect 
on lambda phage replication of several genes associated with 
Fe-S cluster biosynthesis. We performed infection time courses 
of AiscU and other strains deficient in members of the iscRSUA- 
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hscBA-fdx-iscX operon (JiscA, hscB, iscA, iscR, fdx, and iscX). 
The slow growth of the iscS deletion strain made it impossible 
to assess infection dynamics and was not included in the 
analysis. The AiscU, AhscA, and AhscB strains all exhibited 
very similar infection dynamics, clearing the culture signifi- 
cantly before WT (Figure 3 A, top row; Figure 3B). IscU is the 
primary scaffolding protein used for Fe-S cluster biosynthesis 
(Urbina et al, 2001) , and HscA and HscB are accessory proteins 
that help transfer Fe-S clusters to apoproteins (Bonomi et al, 
2008) The AiscA strain exhibited increased replication 
efficiency, but less than the hscA, hscB, and iscUgene deletions 
(Figure 3 A and B) . IscA assists in Fe-S cluster biosynthesis by 
shuttling iron to IscU (Yang et al, 2006) . 

Deletion strains for the remaining three genes in the operon, 
iscR, fdx, and iscX, displayed little if any change in viral 
infection dynamics from WT (Figure 3A and B) . The roles of 
Fdx and IscX in Fe-S cluster biosynthesis are not known 
(Takahashi and Nakamura, 1999), but IscR is the transcrip- 
tional repressor of the ISC operon and its active form includes 
an Fe-S cluster (Schwartz et al, 2001). Interestingly, lambda 
phage replication is slightly but significantly hindered in the 
iscR knockout strain (Figure 3A and B) . Taken together, these 
observations suggest that ISC pathway members IscU, HscA, 
HscB, IscA, and possibly IscR have a more significant role in 
lambda phage infection than Fdx or IscX, and that their WT 
function helps to restrict phage growth. 

We next determined how deletion of tRNA thiolation 
enzymes affected viral susceptibility of the host. We began 
with the Fe-S cluster-containing enzymes encoded by miaB 
and ttcA. MiaB is responsible for methylthio formation of 
ms2i6A at position 37 in tRNA phe (Pierrel et al, 2004); a sulfur 
donor for MiaB has not been identified. TtcA catalyzes the 
2-thiocytidine modification at position 32 of tRNA Arg (ICG), 
tRNA Ars (CCG), tRNA Ar8 (mnm5UCU), and tRNA Ser (CGU) 
(Jager et al, 2004). AttcA strains displayed doubling times 
equivalent to WT E. coli (Figure 3 A, bottom), and little is 
known about the functional consequences of modification. 
Strains lacking either of these genes showed no effect on 
lambda infection dynamics (Figure 3A, bottom; Figure 3C). In 
contrast, deletion strains for the two Fe-S cluster-independent 
tRNA thiolation enzymes Thil and MnmA displayed a 
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Figure 3 The effects of Fe-S cluster and tRNA thiolation deletions on lambda phage infection dynamics. (A) Cell-culture dynamics following infection overlay the 
uninfected control data (WT alone, gray). Small arrows in the cartoon indicate sulfur transfer and the large bent arrow indicates expression of the ISC operon. 

(B) Displacement from the WT strain for ISC operon members (*P<0.05, **P<0.01, bars indicate 95% CI, P-values were calculated using independent two-sample f-test). 

(C) Comparison of the displacement from the WT strain for tRNA thiolation pathway members. Source data is available for this figure in the Supplementary Information. 



significant increase in viral resistance versus the WT strain 
(Figure 3A, bottom; Figure 3C). Sulfur is transferred to Thil 
through direct binding of IscS (Kambampati and Lauhon, 
2000), and Thil then transfers the sulfur to s4U8 of tRNA phe . 
Lambda replication in the Athil strain was mildly but 
significantly inhibited (Figure 3A, bottom; Figure 3C). 

The AmnmA strain exhibited the most dramatic effect on 
viral infection of all the strains associated with tRNA 
thiolation. MnmA catalyzes 2-thiouridination at position 34 



of several tRNAs (Kambampati and Lauhon, 2003) . Sulfur is 
supplied to MnmA through the TUS pathway, beginning with 
the direct binding of IscS by TusA (Ikeuchi et al, 2006). The 
AmnmA strain underwent infection dynamics (Figure 3A, 
bottom; Figure 3C) very similar to the upstream TUS pathway 
gene deletion strains we characterized in a previous study 
(Maynard et al, 2010), including AtusA, AtusBCD, and AtusE. 
All members of the TUS pathway had a strong, significant, and 
similar effect on lambda replication (Figure 3 A and C) , leading 
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Figure 4 Displacement vector comparisons for mutants in £ coli genes known 
to affect frameshifting. Comparison of the displacement from the WT for deletion 
mutants in genes known to affect frameshifting (**P< 0.01 , bars indicate 95% CI, 
P-values were calculated using independent two-sample f-test). Green, reduced 
infectivity; yellow, no effect on infectivity. Source data is available for this figure in 
the Supplementary Information. 

us to conclude that Fe-S cluster-independent tRNA thiolation 
in£. coli exerts a significant effect on lambda phage replication 
through the TUS pathway and, to a lesser degree, Thil. 

Changes in tRNA modification and frameshifting 
propensity inhibit viral replication 

Based on the preceding data and a previous observation of 
increased frameshifting in AmnmA strains (Urbonavicius et al, 
2001), we hypothesized that reduced lambda replication in the 
TUS knockouts was due to increased frameshifting during 
lambda protein synthesis as a result of the decreased tRNA 
modification. To test this hypothesis, we generated infection 
curves for several deletion mutants that have been shown to 
affect frameshifting by loss of tRNA modification, including 
tgt, truA, miaA, thil, and mnmE (Urbonavicius etal, 2001). We 
infected deletion strains for these genes with lambda phage 
and determined their susceptibility to infection (Figure 4; see 
Supplementary Figure 2 for line plots of individual replicates) . 
Atgt, AmiaA, and AtruA exhibited no significant affect on 
lambda phage replication, suggesting that the tRNA modifica- 
tions from these enzymes have no role in lambda phage 
replication. The remaining knockout strains displayed 
significant changes in lambda phage's ability to infect the 
host, including Athil, AmnmA (Figures 3 and 4), and AmnmE 
(Figure 4). Intriguingly, MnmE performs a methylamino- 
methyl modification on carbon 5 of the same tRNA species 
and nucleoside targeted by MnmA (Elseviers et al, 1984), 
suggesting that the antiviral properties of the TUS pathway 
knockouts may be due to specific tRNA modifications instead 
of to general translational fidelity issues. 

Codon usage bias does not underlie differential 
infection dynamics 

Host-preferred codon usage is common among many viruses 
and their hosts (Kunisawa et al, 1998; Lucks et al, 2008). 
We wondered whether the lambda genome used the codons 
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Lysine 


AAA 


0.76 (46116) 


0.65 (486) 




AAG 


0.24 (14174) 


0.35 (262) 


Glutamic acid 


GAA 


0.69 (54431) 


0.57 (515) 




GAG 


0.31 (24 629) 


0.43 (385) 


Glutamine 


CAA 


0.35 (21121) 


0.20 (117) 




CAG 


0.65 (39 836) 


0.80 (459) 



whose tRNAs are modified by the mnmA, and mnmE gene 
products preferentially with respect to E. coli. The tRNA 
species modified by MnmA and MnmE decode the lysine, 
glutamic acid, and glutamine codons (Kambampati and 
Lauhon, 2003); we examined the codon usage for lysine, 
glutamic acid, and glutamine in lambda phage and in E. coli for 
synonymous codon usage bias. Lambda phage encodes lysine 
748 times (65% AAA and 35% AAG), glutamic acid 900 times 
(57% GAA and 43% GAG), and glutamine 576 times (20% 
CAA and 80% CAG), with similar biases in E. coli (Table II). 
This similarity in codon usage does not support the hypothesis 
that lambda phage protein synthesis is any more impaired 
in strains harboring hypomodified 2-thiouridine tRNAs than 
its host. 

Thiolation of tRNA Lys (UUU) and frameshifting via 
PRF are linked through a genetic network 

Since we thought it unlikely that lambda phage replication was 
being significantly affected by general reading frame main- 
tenance in the TUS pathway knockout mutants, we turned 
instead to specific tRNA modifications that could explain our 
observations of an increase in viral resistance. In particular, 
our attention was drawn to PRF as a specific instance of 
frameshifting that could be particularly sensitive to changes in 
frameshifting frequency. 

In lambda phage, PRF occurs during translation of the 
lambda phage GT region, which encodes two possible protein 
products, either gpG or gpGT; expression of gpGT depends on 
ribosomal frameshifting at the slippery sequence GGGTTTG 
(Levin et al, 1993). This sequence encodes the dipeptide 
Gly-Lys in both the 0 open reading frame (ORF; GGT-TTG) and 
the -1 ORF (GGG-TTT). In the 0 ORF, which is maintained 
~96% of the time, gpG is expressed (Figure 5A); and in the 
remaining cases, the ribosomal machinery reaches the slippery 
sequence and slips back one base pair. As a result, the stop 
codon for gpG, immediately downstream of the slippery 
sequence, is bypassed and the larger gpGT product is produced. 
Interestingly, both the lysine codons TTG (0 ORF) and TTT (-1 
ORF) are decoded by the same tRNA Lys (UUU) , a target of MnmA 
and MnmE. We therefore wondered whether TUS pathway gene 
deletion strains would exhibit an altered rate of frameshifting 
and consequently a change in the ratio of gpG to gpGT. 

To monitor gpG and gpGT expression in E. coli, we created a 
plasmid for the inducible synthesis of these proteins (pBAD- 
XGT; Figure 5B) . We PCR amplified the GTregion of the lambda 
genome and introduced this region behind an arabinose- 
inducible promoter, to control expression of the transcript, and 
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Figure 5 Lambda phage proteins gpG and frameshift product gpGT. (A) Schematic of lambda phage's frameshifting region in 6 and T. (B) Schematic of the 
arabinose-inducible region of the pBAD-XGT vector. (C) Immunoblotting of pBAD-M3Tin BW251 13, Aiscil, AtusA, and AtusAAiscU strains. We induced expression of 
the pBAD-A,G7"transcript with 0.2, 0.02, or 0.002% L-arabinose for 2 h and assessed the gpG and gpGT protein levels by immunoblotting against the Xpress Epitope tag. 
We found a decrease in gpGT levels in AiscU and an increase in Aiscil and AtusAAiscU. 



the Xpress Epitope C-terminal tag, to facilitate detection by 
immunoblotting. We then transformed WT, AtusA, AiscU, and 
AtusAAiscU strains with pBAD-XGT. Deletion of the iscU gene 
led to abrogation of frameshifitng, whereas frameshifting 
frequency increased in both the AtusA and the AtusAAiscU 
strains relative to WT (Figure 5C) . We also tested frameshifting 
frequency for additional strains [AmnmA, Athil, AmiaB, 
AttcA, AhscA, AhscB, Afdx, AiscA, and AiscR) and found 
frameshifting levels to be consistent with infection dynamics 
(see Supplementary Figure 3). We interpret these results to 
signify that the decreased thiolation of tRNA Lys (UUU) in the 
TUS pathway knockout strains leads to an increase in 
frameshifting at the slippery sequence. 



Competitive binding of IscU and TusA for IscS 
binding integrates sulfur metabolism and infection 
dynamics 

Surprisingly, the AiscU strain exhibited markedly reduced 
gpGT levels relative to the WT strain (Figure 5C), strongly 
indicating that both IscU and TusA act on lambda phage 
replication by altering the gpG:gpGT ratio. Furthermore, our 
observation that the miaB and ttcA deletion strains exhibited 
normal lambda phage replication (Figure 3) suggested that 
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IscU does not act directly through thiolation of specific tRNA 
modifications. How, then, does IscU influence viral replication 
independent of direct tRNA modification? 

One possible answer is that IscU competes with TusA to bind 
IscS. Recent X-ray crystallography studies showed that IscU 
and TusA bind IscS at distinct but adjacent locations (Shi et al, 
2010). Superposition of the IscU and TusA structures indicated 
a spatial overlap between the volume occupied by TusA and 
IscU when bound to IscS. Furthermore, three-way pull-down 
experiments demonstrated that IscS may bind either TusA or 
IscU individually, but not both at the same time (Shi et al, 
2010). To further explore the possibility that thiolation of 
tRNA Lys (UUU) is affected by competition between TusA and 
IscU for sulfur transfer from IscS, we constructed a mathema- 
tical model of IscU and TusA binding to IscS (Figure 6A) . Our 
model consisted of a set of equilibrium and mass conservation 
relationships. In this case, the model parameters had either 
been previously measured or were not difficult to estimate 
based on existing data (Urbina et al, 2001; Taniguchi et al, 
2010). Our model predicts the relative amount of U34 
thiolation from the TUS pathway for a specified amount 
of IscU. 

We combined this protein binding model with our previous 
mathematical model of E. coli infection by lambda phage 
(Maynard et al, 2010) . Our earlier model predicts the E. coli 
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Figure 6 Competitive binding of IscU and TusA to IscS: theory and experiment. (A) A schematic of our competitive binding model. IscS may bind and transfer sulfur to 
either IscU or TusA, leading to production of Fe-S cluster biosynthesis or thiolated tRNA, respectively. (B-l) Comparison of model predictions (top row) and experimental 
observations (bottom row) of infection time courses. The genetic perturbations include the deletion mutant strains AiscU (B, C) and AtusA (D, E), titration of expression 
from iscil (F, G) and the double deletion mutant strain AtusAAiscU (H, I). (J) The host resistance metric calculated for the AtusA, AiscU, and AtusAAiscU deletion 
mutant strains (**P< 0.01 , bars indicate 95% CI, P-values were calculated using independent two-sample f-test). (K) The prediction of an alternate model— independent 
effect of IscU and TusA on viral biosynthesis— for an infection time course in the AtusAAiscU background. Comparison of this plot with (H) and (I) indicates that the 
competitive binding model is better able to account for the experimental observations. Source data is available for this figure in the Supplementary Information. 



time course dynamics of lambda phage infection based on 
three parameters: the fraction of infections that are lytic 
instead of lysogenic, the rate of infection, and the burst size 
(how many functional progeny are released upon lysis). To 
integrate these two models, we related U34 thiolation (the 
output of the protein binding model) to the burst size (input for 
our infection model) . We were thus able to explore the effect of 
changing the IscU or TusA concentration on the production of 
functional phage and the resulting infection time course. We 
first tested the integrated model's ability to predict the 
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outcome of deleting iscU or tusA, and in both cases the model 
was able to capture the infection time course (Figure 6B-E; see 
Supplementary Figure 4 for line plots of individual replicates) . 
We also considered titration of IscU experimentally as well as 
computationally by introducing an IPTG-inducible iscU gene 
construct into a AiscU strain and determining the infection 
time courses under various IPTG concentrations. The effect 
of IscU titration was predicted well by our competitive 
inhibition model (Figure 6F and G; see Supplementary Figure 
4 for line plots of individual replicates) . 
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Finally, we considered the case of a double deletion of tusA 
and iscU. The competitive inhibition model suggested that the 
double mutant would behave in a manner similar to a AtusA 
strain under infection conditions (Figure 6H). Again, this 
prediction agreed with the experimental time course (Figure 61 
and J; see Supplementary Figure 4 for line plots of individual 
replicates), as well as our experimental observation that the 
ratio of gpGT to gpG in the double mutant was similar to that 
observed in the AtusA strain, but not in the AiscU strain 
(Figure 5C). 

As a control, we compared the predictions of the competitive 
inhibition model with those of the alternative hypothesis that 
IscU and TusA influence viral infection independently. This 
hypothesis required some relatively simple modifications in 
the equations and parameters of our competitive binding 
model (Materials and methods) . Our independent effect model 
was able to capture the results of tusA deletion, iscU deletion, 
and IscU titration as accurately as the competitive inhibition 
model (data not shown). However, the independent effect 
model was not able to predict the behavior of the double 
knockout strain (Figure 6K) . 

Our competitive inhibition model predicts that the removal 
of IscU, one of TusAs primary competitors for sulfur, will 
cause an increase in the fraction of 2-thiolation modified U34 
tRNA Lys (UUU) relative to BW25113. To directly measure the 
fraction of thiolated tRNA Lys (UUU), we performed an 
[(N-Acryloylamino) phenyl] mercuric chloride (APM) north- 
ern blot (Figure 7), which specifically retards the migration of 
the thiolated tRNA fraction (Igloi, 1988). We found a small 
fraction of hypomodified tRNA Lys in BW25113, and as 
predicted by the competitive inhibition model, loss of IscU in 
the AiscU strain decreased this fraction even further. We 
therefore conclude that our competitive inhibition model, but 
not our independent effect model, is sufficient to explain the 
effect of TUS and ISC pathway knockouts on lambda phage 
replication. 



BW25113 AiscU BW25113 AiscU 

P» - tRNA Lys S2U34 

I 



I 

— ■* ' -tRNA Lys U34 

-APM + APM 

Figure 7 APM northern blot of tRNA Lys (UUU) in BW251 13 and AiscU. tRNA 
northern blots for tRNA Lys (UUU) in urea denaturing gels with and without APM. 
tRNA extracts were probed to determine the thiolated fraction of U34 in BW251 13 
and AiscU. A relative shift in the APM + gel reflects the fraction of 2-thiolated 
U34 tRNA Lys (tRNA Lys s2U34). A small but detectable level of hypomodified 
tRNA Lys (tRNA ys U34) was seen in BW251 13 but was essentially undetectable 
in AiscU. 
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Discussion 

Taken together, our observations suggest the existence of a 
complex network that links sulfur metabolism, tRNA mod- 
ification, competitive binding, and PRF-based regulation of 
viral protein expression with host susceptibility to viral 
infection (Figure 8). TusA obtains sulfur from the cysteine 
desulfurase IscS, and passes it along the TUS pathway for 
eventual modification of tRNA Lys (UUU) at U34. The fraction of 
tRNA Lys (UUU) that undergoes this modification influences the 
frequency of PRF in lambda phage due to its role in decoding 
part of the slippery sequence within its GT region. Although 
lambda phage requires both gpG and gpGT for normal 
replication, the gpG:gpGT ratio must be kept high. Deletion 
of any member of the TUS pathway prevents modification of 
tRNA Lys (UUU), and subsequently increases frameshifting, 
both decreasing the gpG:gpGT ratio and lambda phage 
production. 

In the presence of HscA, HscB, and (less critically) IscA, IscU 
binds IscS more strongly than TusA (Shi et cd, 2010) and is 
therefore able to outcompete TusA for sulfur. When IscU, 
HscA, HscB, or IscA is removed from the system, TusA obtains 
more sulfur than normal, leading to hypermodification of 
tRNA Lys (UUU), a decrease in frameshifting, and increases in 
the gpG:gpGT ratio and lambda phage production. Previous 
investigations in B. subtilis reported that disrupting the 
function of IscU affects several levels of tRNA modifications, 
including an increase in 2-thiolation of U34 (Leipuviene et cd, 
2004); our observations suggest that in the case of tRNA Lys 
(UUU), this effect is a result of an increase of sulfur flux 
through the TUS pathway (Figures 5C, 6K, and 7). Removal of 
the transcriptional repressor IscR leads to an increase in 
protein production from the ISC operon (Schwartz et cd, 2001), 
reducing the amount of sulfur available to TusA. This 
competition effect depends on an intact TUS pathway, which 
is why the AtusAAiscU deletion mutant expresses a phenotype 
similar to the AtusA strain. 

This model fits our observations well, but unresolved issues 
remain, most notably the role of Thil in viral infection. Thil 
would seem to be another possible competitor for sulfur from 
IscS (Figure 2A) . However, the Athil strain exhibited decreased 
lambda replication, albeit somewhat less strongly than the 
TUS pathway deletion strains. One clue to this mystery may 
reside with the CyaY protein, which is thought to be an iron 
donor for Fe-S cluster assembly, forming a ternary complex 
with IscU and IscS (Shi et cd, 2010) . The binding site for Thil on 
IscS overlaps substantially with CyaY (Shi et til, 2010) , and it is 
conceivable that in the absence of Thil, CyaY is able to interact 
more efficiently with IscU, either stabilizing the IscSTscU- 
CyaY complex or assisting sulfur transfer to IscU. In either 
case, deleting thil may increase IscU's ability to compete with 
TusA for sulfur. 

Our findings point to several novel antiviral strategies and 
targets. Although this work was performed using E. coli and 
phage lambda, many of the key elements of the network we 
describe are conserved in other systems, in both in host and 
the virus. Both Fe-S clusters and tRNA modifications are 
ubiquitous in cells (Limbach et til, 1994; Py and Barras, 2010) 
and PRF is a common viral translational strategy (Gesteland 
and Atkins, 1996). Indeed, frameshifting has already been 
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Figure 8 A network linking host resistance to viral infection to sulfur metabolism, tRNA modification, PRF, and competitive protein binding. IscS obtains sulfur from L- 
cysteine and passes it to either TusA or IscU. From TusA, the TUS pathway leads to thiolation by MnmA of tRNA Lys (UUU) U34, which is also methylated by MnmE. 
During viral infection, the WT amount of modified tRNA leads to normal translation (and frameshifting) for the viral 6 and 7 genes, and consequently a favorable ratio of 
gpG:gpGT and virion production. Hypomodification of tRNA Lys (UUU), whether by deletion of mnmE or the TUS pathway genes, or by overexpression of IscU, leads to 
increased frameshifting during translation of G and T. This increase in frameshifting leads to more gpGT expression and subsequently, a lower gpG:gpGT ratio and 
decreased virion production. 



identified as a potential target in HIV, but efforts to target the 
gag-pol mRNA secondary structure have been challenging 
(Hung et al, 1998). Our results point to upstream (e.g., tRNA 
modifying enzymes) and even indirect targets (e.g., Fe-S 
cluster biosynthesis pathways), which can just as potently 
alter infection and in some cases are tunable (Figure 6G) . We 
therefore anticipate that our observations and inferences may 
have broad relevance to understanding other infections, 
including viral infection of humans. 

Materials and methods 
Strains 

The E. coli strains used in this study were all taken or derived from the 
Keio Collection, a comprehensive library of single gene deletion strains 
(Baba et al, 2006). BW25113 is the parental strain used for the 
construction of the Keio Collection. Lambda phage was obtained from 
the ATCC (23724-B2). The AtusBCD strain was generated previously 
(Maynard et al, 2010) , and the AtusAAiscU double knockout strain was 
constructed according to the Datsenko-Wanner method (Datsenko and 
Wanner, 2000). Primer sequences are available upon request. The 
lambda GT reporter strain was constructed using a pBAD/His A 
plasmid (Invitrogen, Cat# V430-01). The lambda phage GTregion was 
amplified via PCR (primers: 5'-cacattctcgagatgttcctgaaaaccgaatca-3' 
and 5'-gtgtaaaagcttgccataccggactcctcct-3'). The amplified GTregion 



and pBAD/His A plasmid were cut with Xhol and Hindlll. BW25113, 
AfusA, AiscU, and AtusAAiscU strains were transformed with the 
arabinose-inducible lambda phage GTreporter construct (Figure 5B). 
A chemically inducible AiscU strain was created by transforming AiscU 
with pCA24N-isc(7 isolated from the ASKA collection (Kitagawa et al, 
2005). 



Cell culture quantification 

We monitored infection dynamics using an incubated plate reader 
(Perkin-Elmer Victor3, 2030-0030) under growth conditions described 
previously (Maynard et al, 2010). Exponentially growing cells were 
normalized to 0.1 OD 600 nm. Fifteen microliters of 0.1 OD E. coli and 
15 f-il of ~10 4 plaque forming units (p.f.u.)/ml lambda stock were 
added to 170 yd of LB medium in 96-well plates, representing a 
multiplicity of infection of ~2 x 10~ 4 p.f.u./bacteria. Each plate 
contained four technical replicates of infected and uninfected samples. 
The plate reader measured the absorbance of the culture approxi- 
mately every 17min for 16h. Each plate contained a sample of 
BW25113 for between-plate normalization. Three independent biolo- 
gical replicates were performed for each strain. 



Comparative metrics for time courses 

Time courses were compared based on a useful displacement metric 
that includes a direction in addition to a distance. The distance was 
determined by calculating the Euclidean distance between the two 
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time course vectors; the Euclidean distance between the recorded 
absorbances for each culture were calculated at each time point, and 
these distances were summed for all time points. For our purposes, the 
direction could only hold two values, + /— 1. The positive value was 
associated with strains exhibiting increased resistance to viral 
infection, and the negative value was associated with increased 
susceptibility to viral infection. The direction was determined by the 
difference in the time it took a culture to reach its peak absorbance. For 
example, time courses with longer times to peak compared with 
BW25113 were considered to be displaced in a positive direction from 
BW25113. The displacement was then calculated as [direction) x 
(distance) . 

The Euclidean distance from BW25113 was measured for each 
biological replicate. Statistical significance was determined using a 
two-sample r-test comparing the Euclidian distances calculated for 
each biological replicate and its distance from BW25113 for each plate 
containing the strain being compared. 



E. coli and lambda phage codon usage 

Lambda phage ORFs were extracted from NCBI Reference Sequence 
NC_001416.1 (Enterobacteria phage lambda, complete genome) and 
codon usage was tallied. 

Codon usage for WT E. coli ( W3110) was determined from the Codon 
Usage Database (http://www.kazusa.or.jp/codon/). 



APM northern blot 

The APM northern blot was performed according to the protocol 
specified by Igloi (1988). Briefly, two 8M Urea, 5% acrylamide gels 
were cast in 0.5 x TBE buffer. One gel contained 200 ul of 1 mg/ml 
APM (kindly supplied by Gabor Igloi) for 10 ml of gel. The gels were 
prerunfor 30min at 180 V in 0.5 x TBE. Approximately 2.5 j.ig of RNA 
(for each well) was then mixed with 2 x loading buffer and heated to 
65°C for 3 min before being cooled on ice for 1 min. The gels were run 
for ~ 45 min at 180 V. The RNA was transferred onto a nylon 
membrane using a semi-dry system for 1 h at 250 mA. The transfer 
buffer contained 0.5 x TBE supplemented with 10 mM (3-mercap- 
toethanol. The RNA was then UV crosslinked to the nylon membrane. 
T4 Polynucleotide Kinase (PNK) was used to radiolabel 25 pmol of 32 P 
ATP to 25 pmol of tRNA Lys (UUU) oligonucleotide probe for each blot. 
The oligonucleotide sequence used to probe tRNA Lys (UUU) was 
5'-TGGGTCGTGCAGGATTCGAA-3'. Qiagen's Qiaquick Nucleotide 
Exchange Kit was used to remove free 32 P ATP. The blots were 
incubated in 25 ml of prehybridization buffer (10 x Denhardt's, 6 x 
SSC, 0.1 % SDS) at 42°C for 1 h. The labeled probe was heated to 95°C 
for 2-3 min before being added to the prehybridaztion buffer. The blots 
were hybridized at 42°C overnight. The next day, the blots were 
washed with high salt buffer (5 x SSC, 0.1 % SDS) at 37°C for 20 min 
twice, followed by two washes with low salt buffer (1 x SSC, 0.1% 
SDS) at 42°C for 20 min. After the final wash, the blots were exposed to 
a phosphorscreen for 10 min. 



Assaying lambda PRF 

pBAD-A,GT-transformed cells were cultured overnight at 37°C in LB and 
the appropriate antiobiotics. The cultures were diluted 1:50 into fresh 
LB the following morning and placed back in the shaker for 1 h. The 
culture was then split into 1 ml aliquots arabinose. Arabinose 
induction was performed for 2h at 0.2, 0.02, and 0.002% arabinose, 
after which the samples were spun down in a table-top centrifuge at 
~20 000g for 30 s. The media was aspirated and the pellets were 
frozen at -80°C. The samples were run on a 12% SDS-PAGE gel for 
1.5 h at 80 V in Tris/glycine/SDS buffer. The protein was transferred 
onto a polyvinylidene fluoride membrane for 1 h at 100 V at room 
temperature. The membrane was washed for 30 s in TBST, then 
blocked with 5 % milk blocking buffer for 1 h. The primary antibody 
(Anti-express antibody, Invitrogen) was diluted 1:5000 and incubated 
with the membrane blot overnight at 4 C. The membrane was washed 
five times for 10 min with TBST, followed by a 1-h incubation with a 
1:5000 dilution of the secondary antibody (goat anti-mouse IgG 
horseradish peroxidase). The membrane was washed five times for 
5 min with TBST. Luminol reagent was used for chemiluminescence 
imaging (Santa Cruz Biotechnology, sc2048). 



Mathematical modeling 

We constructed a competitive inhibition model based on the following 
assumed mechanisms: (1) the cysteine desulfurase IscS binds sulfur, 
converting cysteine to alanine in the process; (2) thiolated IscS 
presents sulfur to the apo forms of TusA and IscU; (3) TusA and IscU 
bind thiolated IscS with their respective affinities and obtain the sulfur; 
(4) TusA and IscU transfer the sulfur to their downstream pathways, 
represented as sTUS and sISC, respectively. The source code for our 
model can be downloaded from https://simtk.org/home/lambda-tus. 
The reactions are as follows: 

Cys + IscS «-^-» Cys ■ IscS IscS ■ S + Ala 



IscS ■ S + IscU <-^-» IscS ■ S ■ IscU IscS + IscU ■ S 

IscS ■ S + TusA <-^-> IscS ■ S ■ TusA IscS + TusA ■ S 

fc-s 

IscU ■ S i IscU + sISC 



tRNA enrichment 

BW2511 3 and AiscU were cultured overnight and diluted 1 : 100 to 20 ml 
in the morning. When the cultures reached ~0.5 OD 600nm, the 
samples were placed on ice. The cells were pelleted at 5000r.p.m. for 
5 min at 4°C and the supernatant was decanted. tRNA enrichment was 
performed using standard techniques. Briefly, the cell pellets were 
resuspended in 300 yd cold sodium acetate (0.3 M, pH 4.5, 10 mM 
EDTA) and transferred to eppendorf tubes. The samples were snap 
frozen in liquid nitrogen. After thawing, 300 ul phenolxhloroform (pH 
4.7) was added to the samples and vortexed four times for 30 s bursts 
with 1 min incubation on ice in between. The samples were 
centrifuged for 15 min at 4°C in a table-top centrifuge at maximum 
speed. The top layer (aqueous phase) was transferred to a new tube 
and the acid phenol extraction procedure was repeated. After the 
aqueous phase was transferred to a new tube, the RNA was ethanol 
precipitated by adding three volumes of cold 100% ethanol and 
spinning at maximum speed for 25 min at 4°C. The RNA pellet was 
resuspended in 60 cold sodium acetate (0.3 M, pH 4.5) and 
precipitated again by adding 400 |il 100% ethanol and spinning at 
maximum speed for 25 min at 4°C. The supernatant was decanted and 
all traces of ethanol were removed. The RNA pellet was air dried on ice 
then resuspended in 100 pj cold sodium acetate (10 mM, pH 4.5). 



TusA ■ S -i* TusA + sTUS 

From these reactions we derived the following set of ordinary 
differential equations and mass conservation relationships: 

d[Cys] 



dr 

d[Cys ■ IscS] 



dt 



= fe_i [Cys ■ IscS] - fci [Cys] [IscS] 
fci [Cys] [IscS] - (fe_! + fc 2 )[Cys ■ IscS] 



d[IscS ■ S] 



=fe 2 [Cys ■ IscS] + fc_ 3 [IscS ■ S ■ IscU] + fc_ s [IscS ■ S ■ TusA] 
- fo[IscS ■ S][IscU] - fc 5 [IscS ■ S][TusA] 



d[IscU] 
dt 

d[IscS ■ S ■ IscU 

~~ dt 



fcJlscS ■ S ■ IscU] + MlscU ■ S] - fc 3 [IscS ■ SHlscU] 



= fc 3 [IscS ■ S][IscU] - (fc_ 3 +fc,)[IscS ■ S ■ IscU] 
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d[TusA] 



dt 

d[IscS ■ S ■ TusA 

dT~ 



= fc_ 5 [IscS ■ S ■ TusA] + fc 8 [TusA ■ SI - fc 5 [IscS ■ S][TusA] 



= fc s [IscS ■ S][TusA] - (k_ s + fc 6 )[IscS ■ S ■ TusA] 



d[sISC] 
dt 

d[sTUS] 
dt 



■■ fc 7 [IscU- S] 



= fc«[TusA- S] 



[IscS] + [Cys ■ IscS] + [IscS ■ S] + [IscS ■ S ■ IscU] + [IscS ■ S ■ TusA] 
= [IscS] T0TAL 

[IscU] + [IscU ■ S] + [IscS ■ S ■ IscU] = [IscU] T0TAL 

[TusA] + [TusA ■ S] + [IscS ■ S ■ TusA] = [TusA] total 

[Cys] + [IscS ■ Cys] + [Ala] = [Cys] T0TAL 

Many of the parameter values were taken from the literature. Taniguchi 
et al (2010) reported values of ~ 800 IscS proteins per cell and 450 IscU 
proteins per cell. Approximating the E. coli volume as 10~ 15 L, we 
obtained concentrations of 1.37 and 0.77 uM for IscS and IscU, 
respectively. Since TusA amounts were not reported, we estimated its 
concentration to be similar to but less than IscU 0.5 \iM. An in-vitro 
study reported the if d values of IscS and IscU to be 2 yiM, and the k CM 
[k 2 in the above equations) to be 8.5 per minute (Urbina et al, 2001). 
IscU displaced TusA in a three-way pull-down assay (Shi et al, 2010), 
from which we hypothesized that the K d of TusA-IscU was 20% more 
than that of IscU-IscS. The cysteine concentration was chosen to be 
equal to the IscS concentration. 

The forward rates of reaction were assigned as follows: 
k 1= k 3 =k 5 =k 7 =ks=W 5 , k 4 =k 6 =k_ } . The reverse reaction rates were 
obtained by multiplying the corresponding forward rate by the 
appropriate dissociation constant. Gene deletion conditions were 
simulated by setting the pertinent concentration values to zero. IscU 
titrations were modeled by setting its concentration to 50, 100, and 
150% of the WT amount. 

Steady-state values of the sTUS parameter relative to WT were used 
to scale the amplification factor b in our previous non-dimensionalized 
predator-prey model (Maynard et al, 2010). We simulated infection 
dynamics with the following parameters: fcj=0.25, /=0.93, 
b=1Q . sTUS+sTUSwt ^ j i » =0 .3 ] K*=0A, k s =0. 

The independent effect model hypothesizes that there is no 
competition between IscU and TusA for sulfur. Thus, they each receive 
their own (equally sized) pool of thiolated IscS (iscSu ■ S belongs to 
IscU and iscS A ■ S belongs to TusA) . The dissociation constants of both 
IscU and TusA with their respective pools of IscS were set to 2.7 u.M. 
The concentrations of IscU and TusA were set to the same values as in 
the competitive model. The amount of thiolated IscS, however, was 
increased 100-fold to prevent it from being limiting. The reactions are 
as follows: 



d[IscU] 
dt 



d[TusA] 
di~ 



-■ fc lr [IscU ■ S][IscSu] + fc 3f [IscU ■ S] - ?c lf [IscU] [IscSu ■ S] 



fc 2r [TusA ■ S][IscS A ] + fc 4f [TusA ■ S] - fc 2f [TusA][IscS A ■ S] 



d[IscSu ■ S] 



At 



d[IscS A ■ S] 
dt 



= fc lr [IscU ■ S][IscSu] - fc lf [IscU] [IscSu ■ S] 



= fc 2r [TusA ■ S][IscS A ] - fc 2f [TusA] [IscS A ■ S] 



d[sISC] 
dt 



d[sTUS] 



dt 



fc 3f [IscU ■ S] 



: fc, f [TusA- S] 



[IscU] + [IscU ■ S] = [IscU]. 



[TusA] + [TusA ■ S] = [TusA], 



[IscSu ■ S] + [IscSu] = [IscSu ■ S]. 



[IscS A ■ S] + [IscS A ] = [IscS A ■ S] T 



. • = fIscS A ■ SI- 

The kinetic parameters were specified as follows: 
fc lf =fc 2f =fc 3f =fc 4f =10 5 . Non-zero reverse rates were set to the product 
of the forward rate and the appropriate dissociation constant. 

Steady-state values of the sTUS and sISC parameters relative to 
WT were used to scale the amplification factor b. Parameters 
for the predator-prey equation were ^=0.25, f=0.93, 
b=10 ■ sT " s T nc TUSwT - 50 ■ slsc ,^' SCwT , li * — 0 . 3 , K* =0.4, fc s =0. We see that 
in both models, the WT burst rate is 20 and the tusA knockout burst 
rate is 10. Supplementary Table 1 describes the parameters used in the 
two models. 

The standard ode23s solver in Matlab was used for solving the 
sulfur transfer equations. The ode45 solver was used to solve the 
predator-prey equations. 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 



IscU + IscSu ■ S IscU ■ S + IscSu 



TusA + IscS A ■ S TusA ■ S + IscS A 



IscU ■ S -=> IscU + sISC 



TusA ■ S -!=> TusA + sTUS 

From these, we derived the following set of ordinary differential 
equations and constraints: 
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